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Abstract 


Hexagonal boron nitride (hBN) is a natural hyperbolic material which can also accommo¬ 
date highly dispersive surface phonon-polariton modes. In this paper, we examine theoreti¬ 
cally the mid-infrared optical properties of graphene-hBN heterostructures derived from their 
coupled plasmon-phonon modes. We found that the graphene plasmon couples differently 
with the phonons of the two Reststrahlen bands, owing to then - different hyperbolicity. This 
also leads to distinctively different interaction between an external quantum emitter and the 
plasmon-phonon modes in the two bands, leading to substantial modification of its spectrum. 
The coupling to graphene plasmons allows for additional gate tunability in the Purcell factor, 
and narrow dips in its emission spectra. 


1 Introduction 

Polaritons are hybrid modes of photons and charge dipole carrying excitations in crystals. Two of 
the most common types are surface plasmon polaritons 1-3 and phonon polaritons. 4-7 Such modes 
have been shown to be of technological relevance in subwavelength imaging, 8,9 biosensing, 10,11 
waveguiding, 2,12 photovoltaics 13 and quantum information. 14 16 While surface plasmons rely on 
free electron oscillations, surface phonons exist because of the lattice vibrations in polar crystals. 

Graphene has been shown to be good candidate for tunable plasmonics in the mid-IR and tera¬ 
hertz range, 16-21 owing to the possibility of electrostatic doping 22 and its ability to produce higher 
confinement and lower losses compared to metals. 23 On the other hand, recently near-field imag¬ 
ing has shown that phonon polaritons in hexagonal boron nitide (hBN) possess extremely high 
confinement and even lower loss compared to graphene plasmon polaritons. 5 hBN shows natu¬ 
ral hyperbolicity 7 which can potentially be used to explore exotic photonic properties 24 such as 
strong spontaneous emission enhancement, 25,26 negative refraction 27 and thermal radiation en¬ 
hancement. 28 Since both graphene plasmons and hBN phonons reside in the mid-IR, the optical 
properties of graphene-hBN heterostructures would allow one to marry the advantage of their con¬ 
stituents, electrical tunability in the former and high quality factor of the latter, through their hybrid 
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plasmon-phonon polaritons. 

The study of the optical properties of graphene-hBN heterostructures are also motivated by the 
following recent developments. Firstly, hBN is now being used as a substrate of choice due to 
the preservation of high carrier mobility, as opposed to conventional SiC>2 substrates. 29 The higher 
carrier mobility also translates to better plasmon quality factors. 30 Secondly, phonon modes of 
hBN can couple to graphene plasmon providing the possibility of observing interesting effects 
such as phonon-induced transparency. 31,32 Recently, a study of patterned graphene on monolayer 
hBN revealed a coherent coupling between plasmon modes in graphene and optical phonon modes 
in single layer hBN. 33 As shown in, 5 hBN thin films can support several higher order phonon- 
polaritonic waveguide modes inside the Reststrahlen band. These modes show a dispersion which 
can be efficiently controlled by varying the thickness of the slab. In addition, both graphene and 
hBN can be grown in large area using CVD techniques, 34,35 hence allowing in principle the con¬ 
struction of arbitrary heterostructures multilayers stack. The combination of these properties ren¬ 
der graphene-hBN heterostructure an interesting photonic system. 

£ Graphene 

* / 



(a) (b) 

Figure 1: (a )Schematic of the geometry. Monolayer graphene deposited on an hBN film of thickness thBN- The figure is 
not to scale. This work focuses on the interaction of graphene plasmon with phonon in hBN (left) and the emission of 
a nearby dipole into these hybrid modes (right). (b)Permittivity of liBN: Permittivity tensor components of hBN clearly 
show the possibility of hyperbolicity in certain frequency ranges. The relevant parameters for the phonon frequencies 
were obtained from. 36 


Typical configuration of the graphene-hBN heterostructure studied in this work, that is, mono- 
layer graphene deposited on hBN thin film, is shown in Fig. la. First we will describe the nature 
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of the couplings between graphene plasmons and phonons of hBN for two different regimes of the 
hyperbolicity of hBN. Secondly, we exploit this coupling to demonstrate the possibility of inducing 
a dip in the Purcell spectrum due to this plasmon-phonon coupling and show its tunability via the 
external knobs of hBN slab thickness and of active tuning using electrostatic gating. 


2 Optical response of hBN and graphene 


Hexagonal boron nitride: Hyperbolic materials are anisotropic materials where the relative per¬ 
mittivity tensor is such that one of the three components has a sign different from the other two. 
This property leads to a hyperbolic or indefinite dispersion for electromagnetic waves propagating 
inside such a material, which results in exotic photonic properties. 24 Until recently, most practical 
realizations of hyperbolic media relied on artifically engineered systems, or the so called metama¬ 
terials. However, the recent discovery of a natural hyperbolicity in hexagonal boron nitride (hBN) 
crystals would make possible the design of atomic scale hyperbolic metamaterials, and might po¬ 
tentially allow one to cross over to regimes beyond the simple effective medium description in 
conventional hyperbolic metamaterials. 

hBN is a van der Waals crystal with two kinds of IR active phonon modes relevant to hyper¬ 
bolicity: 1) out of plane A 2 U phonon modes which have cojo = 780cm -1 , a>io — 830cm” 1 and 
2) in-plane E\ u phonon modes which have g>to = 1370cm , CQlo = 1610cm” 1 . 36 This leads to 

two distinct Reststrahlen (RS) bands, where the lower frequency RS band corresponds to type-I 
hyperbolicity (£|| < 0. £± > 0) and the upper RS band shows type-II hyperbolicity (£_ < 0, £|| > 0). 
The hBN permittivity is given by 


— (t Oto, mY 

(®TO,m) 2 ~ (O 2 — l(OT m 


( 1 ) 


where m =_L, ||. The parameters employed in the above equation are taken from. 36 In addition 
to the LO and TO frequencies which have been mentioned in the previous paragraph, the other 
parameters are = 4.87, £«, n = 2.95, r_|_ = 5 cm” 1 and F|| = 4 cm” 1 . 
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Graphene: The graphene response is modeled using local random phase approximation (local 


RPA). At temperature T, the two dimensional conductivity of graphene is given by: 
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where H(a>,T ) = smh(fKo/kT)/[cosh(p/kT) +cosh(ha>/kT)}. The first term in Eq. (??) repre¬ 
sents intraband contribution and the remaining terms are contributions of the interband transitions 
to the total graphene conductivity. Here T is the electron relaxation time. While Landau damping 
itself is already included in the conductivity model, the relaxation time typically has other con¬ 
tributions from 1) impurity scattering, 2) scattering with phonons (tiCQoph = 0.2 eV) in graphene 
and phonon modes of polar substrates, 3) higher order processes such as phonon coupled to e-h 
excitations (which have to be treated separately), etc. 20,30,38 In literature, relaxation times as long 
as 1000 fs have been reported. 29,39 For frequencies larger than the optical phonon frequency of 
graphene, typically r ~ 50 fs. 20 In this work, we use a graphene DC mobility p D c of 10,000 
cm 2 V 's-'. The mobility and relaxation time are related by z = pix'E p / evj . 23 


3 Plasmon-phonons: two regimes of mode coupling 

We consider the geometry as shown in Fig. la. In this general case, by invoking the quasistatic 
approximation, the modal dispersion inside the two RS bands can be written as: 


q{(0) 
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thBN 
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( 3 ) 


where y/ = ^£||/£_l/z, o is the conductivity of the graphene, is the thickness of the hBN film, 
£ a and £ s are the relative permittivities of air and the substrate. On the other hand, outside the RS 
bands, the above condition is not sufficient since t/r becomes imaginary. Thus, outside the two RS 
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bands, the modal dispersion takes the form: 


1 ln W ^-iY{e±/£g)+i(q/k 0 ){Z 0 a/£ a ) \ _ / 1 -iy(£±/£ s ) \ 

2 thBN n [ V 1 + iv(£±/£ a ) +i(q/k o) (z 0 <7 /£ a )) \ 1 + W(^/^) ) 


( 4 ) 


In the lossless case, inside the two RS bands, y is real whereas ty > 0 outside the RS bands. 
It should be noted that for all subsequent calculations, we have used full-wave calculations, but 


we will refer to the quasistatic dispersion equations given above for intuitive understanding (see 
supplemental information for a comparison between quasistatic and full-wave results). 




hBN (AHA) 


Graphene-hBN (AGHA) 


Figure 2: Dispersion curves and mode profiles for phonon-polaritons in hBN thin slab (AHA) and coupled phonon- 
polariton plasmon modes in graphene-hBN system (AGHA): Subfigures (a), (b), (e) and (f) represent the dispersion 
and mode profiles near the in-plane resonance; (c), (d), (g) and (h) represent those near the out-of-plane resonance 
of hBN. The alphabetical labels at different (k x . (o) points on each of the dispersion diagrams are independent of 
each other and only correspond to the field profiles on their immediate right. Thickness of hBN is tknN = 50 nm and 
graphene Fermi level is assumed to be 0.5 eV. Here AHA and AGHA denote air-hBN-air and air-graphene-hBN-air 
respectively. 


To understand the plasmon-phonon coupling, we consider two cases: a) air-hBN-air (AHA) and 
b) air- graphene- hBN-air (AGHA). A third case of air- graphene- hBN- graphene- air (AGHGA) 
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has also been presented in the supplementary information. We have chosen the geometry of sym¬ 
metric waveguides (the two surrounding half-spaces are both assumed to be air) only for simpler 
presentation, but all our results and discussions in this section have been verified and hold true for 
asymmetric waveguides, which is the common situation experimentally. Later in this paper, when 
we do consider a substrate, these geometries will be denoted as AHS, AGHS and AGHGS instead. 
The refractive index of this dummy substrate was chosen to be 1.5. 

The phonon-plasmon-polariton dispersion is highlighted in the frequency range close to the 
two RS bands in Fig. 2e and 2g. The nature of these modes is completely different between the 
in-plane and out-of-plane RS bands owing to differing types of hyperbolicity. This is described 
next two subsections. Different values of k x on these dispersion diagrams can be experimentally 
accessed for instance by fabricating finite width graphene nanoribbons on top of hBN. 20 All the 
dispersion curves are represented by a density plot of the imaginary part of the Fresnel reflection 
coefficient for p-polarization in (co,k x ) space. The field profiles were obtained using COMSOL 
MULTIPHYSICS package. 

3.1 Plasmon-phonons for the in-plane RS bands 

This is the frequency range which shows type-II hyperbolicity, which means that the in-plane 
dielectric function is negative and the out of plane component is positive. Therefore, for this band 
the quantity y in Eq. ?? is less than zero also. For the AHA case, this gives several slab polariton 
modes for n = 0,1,2,3,... as suggested in the dispersion in Fig. 2a. Each of these mode numbers 
corresponds to the number of nodes in the tangential electric field E x as shown in Fig. 2b. It 
is apparent from Fig. 2b that the fields are mostly confined inside the hBN slab rather than the 
surrounding media or the interfaces. It should be noted that the modes with smaller number of 
nodes occur at higher frequencies. This is a consequence of type II hyperbolicity and the shape of 
the permittivity curve for the hBN slab. 

When the top interface of the hBN slab is covered with graphene (AGHA case), we observe 
an interesting effect: the n = 0 phonon polariton merges smoothly into a plasmon mode of the 
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graphene as shown in Fig. 2e. Note that the phonon polariton modes inside the hBN are sinusoidal 
in the z direction, whereas the plasmon mode is exponentially decaying. This is also apparent from 
a comparison of the field plot (Fig. If) for frequency points labeled “a" and “c" between the AHA 
and the AGHA geometries where in the latter case, the field profile inside the hBN is deviating 
markedly from a sinusoidal profile and tending toward an exponentially decaying one. This has 
the additional effect of shifting the peak of the field from the bulk of hBN to the graphene. 

This smooth transition is unexpected especially given the fact that at the top of this RS band, 
the hBN undergoes a topological transition 40 from hyperbolic to elliptical dispersion. However, 
an abrupt change at such a transition is expected only in bulk crystal of hBN that is thick enough. 
Since we have a thin layer of hBN, this transition can be smoothened. This behaviour can also 
be understood in terms of a mode coupling picture. The graphene plasmon dispersion below this 
RS band starts becoming phonon-like as it approaches the TO frequency of the hBN. On the other 
hand, the lowest (zero) order phonon-polariton mode acquires a plasmonic character near the LO 
frequency. Higher order phonon polariton modes are not affected by this coupling due to symmetry 
mismatch which can be seen by looking at the field profiles of the modes “d" and “e" in Fig. 2f. 
The plasmon mode “f" does not show any node inside the hBN whereas all the phonon-polariton 
modes other than “c" show one node or more. Thus mode “c" which is the lowest order mode can 
transform smoothly into the plasmon. This hand waving argument can be rigorously substantiated 
by calculating the modal overlaps. 

In order to confirm this picture, we investigated whether by coupling to the acoustic plasmon 
mode of the double layer graphene we can make the first order phonon-polariton mode “d" also 
acquire plasmonic character. This case (AGHGA) is presented in the supplementary information, 
where it is clearly seen that the first order phonon polariton also merges smoothly into the acoustic 
plasmon mode. This behaviour can be easily understood by noticing that the graphene plasmon 
has split into symmetric and antisymmetric modes with zero and one nodes which will couple only 
to the zeroth and first order phonon polariton modes respectively. 
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3.2 Plasmon-phonons for the out of plane RS bands 


Phonon polariton modes in the out of plane RS bands for thick hBN slabs have received relatively 
less attention in literature so far since the out of plane phonon polariton response is relatively 
weak. 41 This band displays type I hyperbolicity where the out of plane dielectric function is neg¬ 
ative instead of the in-plane one. This difference makes its coupling to the graphene plasmon 
markedly different from the in-plane phonon case discussed above. 

For this band, y/ in Eq. ?? is positive. Thus, in order to make the real part of the wavevector 
q positive, we need n < 0. This is also reflected in the field plots in Fig. 2d where it is observed 
that there is no mode with zero number of nodes. In other words, the minimum value of \n\ is 

I instead of zero. A simple explanation of this behaviour can be obtained by reference to the 

relative signs of the Poynting vector and the wave-vector. 42 In the case of out of plane resonance, 
since the dispersion is type-I, the Poynting vector component S- and the wavevector component k z 
point in the same direction. Thus, the dispersion equation k z thBN + does not admit n — 0 

solution since the left hand side is strictly positive. On the other hand, for the in-plane mode, type- 

II dispersion causes S z and k : to point in opposite directions which permits the left side approach 
zero value, thus allowing n = 0 solution. Such behaviour also consistent with the predictions in. 42 

In this case one can observe that the ordering of the modes in the region with negative group 
velocity is such that that the lower order mode occurs at smaller frequency, unlike the case of 
in-plane RS band. Again, this is due to the difference in the type of hyperbolicity. 

In the AGHA case, it is observed that below the RS band, the plasmon mode becomes phonon 
like as it approaches the TO phonon frequency. The plasmon branch on the top of the RS band ac¬ 
quires a phonon-polariton character near (Qlo- We also observe the effect of mode repulsion which 
causes the first order phonon polariton to blueshift compared to the case without graphene. As 
shown in the supplementary, this behaviour persists for the AGHGA case where both the graphene 
plasmon modes acquire phonon-polariton like character near ft)/ o and (Ojo- 

In summary, we observe that the coupling between the graphene plasmon and the hBN phonon- 


9 


polaritons is completely different between the two RS bands. This is attributable to the different 
types of hyperbolicity associated with the in-plane and the out-of-plane phonon modes. In the next 
section, we will attempt to use this coupling for tuning the local density of states. 

4 Tunable spectral dips in emission 

When an excited quantum emitter is placed near a material system which can support a photonic 
mode, its lifetime is modified compared to the case when it emits in free space. This phenomenon 
is called Purcell effect. 41 The emitter can typically be an excited atom, molecule or quantum dot. 

The emitter can release its excitation energy into free space radiative modes as well as resonant 
modes and non-radiative or lossy modes of the material. The radiative contribution is modified due 
to the change of the boundary conditions for the electromagnetic fields because of the presence of 
the material. This results in the modification of far-field emission patterns. 44 

More interestingly, the emitter can also release its energy through the available resonant modes 
in the neighbouring matter. The strength of this light-matter interaction is governed by the ratio of 
quality factor and the volume of the available photonic mode. The basic idea of nanophotonics is to 
provide sub-wavelength mode volumes which enhance this interaction. For instance, the power of 
graphene plasmonics is in providing very small mode volumes. 16 On the other hand, hBN phonon 
polaritons provide in addition, high quality factors as well. 5 

One key aspect of hyperbolic materials is that when moving from a hyperbolic dispersion 
regime to an elliptical one, a sharp change in the local density of states (LDOS) is observed. 40 
This is due to the availability of high-k states in the hyperbolic band. To explore this phenomenon 
in the context of graphene-hBN slab system, we considered the spontaneous emission enhancement 
(Purcell effect) of a quantum emitter with polarization perpendicular (z.) to the hBN surface. This 
orientation of the emitter was chosen since the z-dipolc moment will only couple to the TM modes. 
Parallel polarization can also be treated in a similar manner. In the present case, the spontaneous 
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emission rate of the z-polarized dipole, also called the partial LDOS or PLDOS, is given by: 40,45 


r 

To 


l + 2k o- 


dk T k x 



(5) 


where To is the free space radiative decay rate, d s is the distance of the quantum emitter (source) 
from the hBN slab and r p ((0,k x ) is the Fresnel reflection coefficient from the graphene coated hBN 
slab. This expression includes both radiative as well as non-radiative contributions to the total de¬ 
cay rate. The Purcell enhancement in the hyperbolic regime is clearly observable in the AHS case 
shown in Fig. 3a and 3b, where the PLDOS rises sharply inside both the RS bands. As shown in 
the following, in the presence of graphene, a dip in the Purcell enhancement is observed, whose 
spectral width and location is strongly tunable both actively using electrostatic doping and also by 
changing the thickness of hBN. In particular, for the out-of-plane RS band, the reduction in the de¬ 
cay rate is consistently found to be about an order of magnitude in the presence of phonon-plasmon 
coupling. The observed dips in Purcell spectrum arise are due to coupling between spectrally broad 
plasmon and much narrower phonon resonance. Such dips is analogous to induced transparency in 
the absorption spectrum observed recently in such systems. 46 The induced transparency can be un¬ 
derstood in the classical coupled oscillator picture with oscillators of contrasting damping rate. 32,46 
Such a large modification in the PLDOS suggests possible application in tuning the radiative en¬ 
ergy transfer to hBN via phonon modes. Note that the presence of graphene is also expected to 
enhance the nonradiative decay rate of the emitter via energy transfer into the plasmon mode. 16 


4.1 Purcell spectra versus emitter position 

The shape of the Purcell spectra in this case strongly depends on the dipole distance as can be 
seen in Fig. 3a and 3b. To understand this variation we note that since we are in the near field, 
the integrand in Eq. ?? can be simplified as °c k x e~ 2/i ’ y/ '3 { r p (o). k x )}. As noted in, 47 this form 
of the integrand leads to the maximum contribution to the PLDOS coming from k x « 1 jd s . This 
effective wave-vector cutoff together with the dispersion curve can be used to predict the shape 
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Wavenumber cm 

(a) AHS 



Wavenumber cm 


(b) AHS 




Wavenumber cm 1 Wavenumber cm 

(c) AGHS ( d ) AGHS 


Figure 3: Dipole distance dependent Purcell spectra in graphene-hBN system: Distance d s of the quantum emitter 
from the hBN is varied. Subfigures (a) and (b) represent the spectra for the case of hBN thin film (AHS) for the two 
RS bands respectively. Subfigures (c) and (d) represent the case of monolayer graphene deposited on the hBN thin 
film (AGHS). Dashed lines show the Purcell enhancement due to graphene plasmon without the phonon-polariton 
contribution from the hBN slab, that is, assuming —>• 0). Thickness of hBN slab is fixed at f/,&y = 50 nm. 

Graphene Fermi level is assumed to be 0.5 eV. In all the figures, the emitter is on the air side. 
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of the Purcell spectra at different distances of the dipole emitter. For instance, if we look at the 
spectra near the out-of-plane modes in Fig. 3a, it is seen that as the dipole moves away from the 
hBN, the PLDOS peak shifts towards the upper part of the RS band, that is, towards (Qlo- To 
understand this, we note from the dispersion curve shown in Fig. 2c that at high k x , the dispersion 
of the phonon-polariton modes becomes flatter, implying a high density of states available at the 
respective temporal frequencies. Thus at small d s , all of these states are contributing to the PLDOS 
which remains high throughout the RS band. As the emitter is taken farther away from the the 
hBN, the maximum contributing k x moves from higher to lower values. At smaller values of k x , 
the lower phonon-polariton branches start acquiring curvature which results in reduction of density 
states. Thus at these values of k x (or equivalently, d s ) only the high frequency phonon-polaritons 
which have a flatter dispersion, can contribute effectively to the PLDOS. This results the in PLDOS 
peak becoming sharper and moving upward in frequency as the dipole moves away from the hBN. 
A similar behaviour is observed for the in-plane phonon RS band in Fig. 3b, where due to the 
opposite shape of the dispersion curve, the PLDOS peaks move downward in frequency as the 
dipole separation from hBN increases. 

When graphene is coated on top of this hBN slab, even more interesting spectral shapes are ob¬ 
tained. Most notably, due to the plasmon-phonon coupling, spectral dips of an order of magnitude 
can be obtained in the PLDOS for a range of dipole distances. Let us first focus on the out-of-plane 
RS band for the AGHS case, as shown in Fig. 3c. It is clear from the dispersion in Fig. 2g that 
at large k x , the plasmon dispersion below (Ojo is flat and this results in a peak in the PLDOS for 
smaller d s in Fig. 3c. This plasmon peak peak seems to disappear as we move to larger d s . This can 
be explained based on two observations. Firstly, since the plasmon dispersion acquires curvature 
at low k x , there is a broadening of the plasmon peak since it is no longer flat. Secondly, at low k x , 
the plasmon branch in the dispersion in Fig. 2g redshifts. Both these effects lead to an apparent 
disappearance of the PLDOS peak near (Ojo with increasing d s . 

Next, let us try to understand what happens when we enter the RS band from below 0)/ o- If we 
are at sufficiently large d s (or low k x ), the plasmon is far below cojo an d it has curvature, implying 


13 


that it will contribute a broad peak located far below (Ojo■ Inside the RS band on the other hand, 
we know from previous analysis that there are flatter states available. Thus these two effects give 
the impression of a kind of “band gap" being created near (Qrch which will result in a decrease in 
the PLDOS in between these two peaks since there are no states available. Note that the spectral 
width and location of this dip will depend on the value of d s that governs which k x contributes 
maximally. As shown in the supplementary, we can even make this “dip" go away depending on 
the dipole distance d s . The width of this dip is smallest at large k x (small d s ), which is explained 
by the changing curvature of the phonon-polariton branches as a function of k x , as explained in the 
AHS case before. This explanation is consistent with our findings in Fig. 3c. 

Beyond (Qlo, we encounter the plasmon mode again. Unlike the plasmon mode below (D/ q, 
this plasmon mode is flatter at low k x because it acquires phononic character there. Thus the 
plasmon peak above the RS band is narrowest for large d s . This explains the trend in the PLDOS 
peak arising above CO 10 in Fig. 3c. 

This spectral signature is carried over in a qualitatively similar way to the in-plane RS band as 
shown in Fig. 3d with notable differences in the trend peak locations on account of the different 
type of hyperbolicity of this band. For instance, the dip below (Oto in this case is always outside 
the RS band. This is because the higher order phonon-polariton modes for the in-plane branch 
cluster near (Ojo instead of (Olo like in the out of plane case. This results in there being always 
states available at frequencies even slightly above COjo- 

4.2 Purcell spectra versus hBN thickness 

The dispersion of guided modes in the hBN slab depends strongly on its thickness. This variation 
in the dispersion has a profound influence on the PLDOS spectra. Let’s focus on the AHS out of 
plane phonon case first. In this case, as shown in Fig. 4a, with increase in hBN thickness 
the phonon-polariton dispersion moves closer to (Ojo or the lower end of the RS band at small k x . 
The implication of this movement is that the PLDOS peak starts moving from the upper end of the 
RS band (at small thBN), broadening toward the lower end (at large Ihun)- The broadening occurs 
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(a) Dispersion trends: out of plane 



Wavenumber cm' 


(b) AHS 



Wavenumber cm 


(c) AGHS 



(d) Dispersion trends: in plane 



Wavenumber cm 


(e) AHS 



Wavenumber cm 


(f) AGHS 


Figure 4: hBN thickness dependent Purcell spectra in graphene-hBN system: (a) and (d) represent trends for dispersion 
variation as a function of hBN thickness and graphene doping near the out of plane and in-plane phonon resonances. 
We consider here the AGHS system. The arrows point in the direction of increasing thm or Ep. The trends shown 
are for independent variation of the two variables and not simultaneous. Unlike other dispersion plots in this paper, 
here the v-axis represents k x to make it easier for the reader to relate it to 1 /d s . See text for details, (b), (c), (e) and 
(f) represent hBN thickness dependent Purcell spectra in graphene-hBN system. Distance d s of the quantum emitter 
from the hBN is fixed at lOOnm. In the AGHS case, the emitter is on the graphene side. Dashed lines show the Purcell 
enhancement due to graphene plasmon without the phonon-polariton contribution from the hBN slab, that is, assuming 
► 0). Graphene Fermi level is assumed to be 0.5eV. 
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again because as t hBN is increased, 1) more states are available throughout the RS band (at fixed k x ) 
and 2) the redshifted phonon polaritons, specially at low k x , have greater curvature. Both of these 
expectations are verified in the calculations shown in Fig. 4b. 

With the graphene coating on the hBN slab, there is a spectral dip immediately below COpo- As 
explained for the d s variation earlier, this is due to the broad peak arising from the curved nature 
of graphene dispersion at the chosen d s , which is peaked much below copo- As t hBN increases, the 
width of the dip becomes narrow because it is “pushed" from the right by the phonon-polariton 
states which are being red-shifted. This explains the movement of the spectral dip in Fig. 4c. 

This behaviour persists for the in-plane modes as shown in Figs. 4e and 4f, with the difference 
being the opposite trend of peak movements, on account of differing types of hyperbolicity, as 
explained earlier. 


4.3 Active control of spectral dips 




Wavenumber cm 1 Wavenumber cm -1 

(a) AGHS (b) AGHS 

Figure 5: Fermi level dependent Purcell spectra in graphene-hBN system: The doping of the graphene is varied. The 
thickness of the hBN slab is 50 nm. Distance d s of the quantum emitter from the hBN is fixed at lOOnm. The AGHS 
case has been considered with the emitter is on the air side. Subfigure (a) represents Purcell spectra around out-of- 
plane resonance and (b) represents that near in-plane resonance. Dashed lines show the Purcell enhancement due to 
graphene plasmon without the phonon-polariton contribution from the hBN slab, that is, assuming e/,^(ft) —» 0). The 
inset in subfigure (b) presents the decay rate near co = 1610cm 1 as a function of graphene doping. 

The Fermi level of graphene can be tuned via electrostatic doping. 48 This is one of the main 
motivations for using graphene for photonics applications. 18,49 Importantly, graphene with higher 
doping has a plasmon dispersion which is closer to the light line. 16,17,50 Thus for a fixed dipole 
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distance and hBN thickness, changing the Fermi level offers a route for active control of the spectral 
location and width of the spectral dips in PLDOS. 

To demonstrate such a gate tunability of the PLDOS, we carried out calculations for the AGHS 
case as shown in Fig. 5. For the out of plane case, we expect that at larger doping (or larger Ep), 
the graphene plasmon below (Oto will flatten out at smaller k x , as shown in Fig. 4a. Conversely, 
at a fixed value of d s , this graphene plasmon branch will have more curvature and will be farther 
below (Oto f° r smaller doping. On the other hand, due to availability of states inside the RS band, 
the PLDOS is expected to rise there. This results in a dip occurring between the two peaks near 
(Oto- 

Secondly, towards the upper end of this (out of plane) RS band, we observe an interesting 
trend. The half line width of the part of PLDOS peak which is inside the RS band decreases with 
increasing Ep. On the other hand, the line width of the PLDOS peak outside the RS band (at the 
higher frequency end), increases with increasing Ep. There are two effects in play here. Firstly, for 
smaller Fermi-level, since the graphene plasmon is very highly confined to the graphene itself, it 
does not interact very strongly with the phonon-polariton. As the doping is increased, the plasmon 
confinement decreases and as shown in Fig. 4a, the phonon-polaritons inside the RS are repelled. 
At reasonably small k x , not only does this repulsion cause a blue shift of the phonon-polaritons 
but also reduces the curvature, resulting in narrower linewidth inside the RS band at high graphene 
doping. For frequencies above COpo , the graphene plasmon is flatter for smaller doping. This causes 
the linewidths of the PLDOS peaks to the right of (Opo to broaden as doping is increased. A subtle 
point to note here is that the graphene damping is expected to decrease with increasing Ep since 
we assumed the relaxation time z to be proportional to Ep. But we have also checked that these 
trends persist for the case of constant z. Similar but complementary trends are observed for the 
in-plane modes as shown in Fig. 5b. Note that due to high frequency, the interband excitations 
begin to damp the graphene plasmon mode, which results in an overall decrease in the PLDOS of 
the uncoupled graphene itself. Interestingly, for the in-plane band, the spontaneous emission rate 
is most sensitive to graphene doping near (Opo (1610cm -1 ). On the other hand, further inside this 
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RS band at lower frequencies, the decay rate is almost independent of the graphene Fermi level. 
Therefore, to study the emission near ft)/o, in the inset of Fig. 5b, we have plotted the decay rate 
at (Qlo , as a function of the graphene Fermi level. First of all, we observe that as we go from 
higher to lower doping, the decay rate becomes smaller. This trend can be explained as before, 
using Fig. 4d. In accordance with Fig. 4d, with decreasing Ep, the phonon polariton plasmon 
modes redshift, resulting in lowering of the decay rate at this frequency (note that we have to look 
at constant k x — I fd s ). Conversely, one can say that for this emitter location, higher gate voltage 
gives higher decay rate which arises due to a larger contribution from the phonon polariton modes. 

However, the above picture explains the trend only upto a certain doping. As shown in the 
inset of Fig. 5b, below a certain Fermi level, the decay rate starts rising again and achieves a con¬ 
stant value corresponding to undoped graphene. This critical value occurs near Ep ~ ft), which 
in this case happens to be around Ep — 0.2eV. To explain this, we note that as the frequency 
approaches the Fermi level, while the plasmon mode becomes weaker, the interband contribution 
becomes larger. We have checked that if we remove the interband damping term from the graphene 
conductivity, the decay rate continues to decrease as the Fermi level becomes smaller (see supple¬ 
mentary information). This suggests that the rise of the decay rate at small doping is explained by 
a decay route opening up through this lossy channel of interband transitions as we asymptotically 
approach the case of undoped graphene. 51,52 At zero temperature, this critical point should occur 
near Ep ps ffl/2. However due to finite temperature smearing in the graphene conductivity, it is 
observed to occur near Ep « ft). On decreasing the temperature we observed that the critical Fermi 
level indeed approaches ft)/2, as shown in the supplemental information. 

Lastly, it should be noted that in this case, the shape of the decay rate versus Ep (as shown in 
the inset of Fig. 5b) is dependent on the emitter location (d s ). In general we expect there to be an 
optimum value of Ep where the decay rate is maximized, corresponding to the plasmon-phonon 
dispersion coinciding at the given frequency ft) and wavevector k x ~ 1 /d s . In the supplementary 
information, we further study the shape of this curve for different emitter locations. 
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5 Conclusion and outlook 


In this work, we have shown that heterostructures comprised of graphene and hBN can provide a 
highly versatile tool for controlling light matter interaction at the nanoscale. The gate tunability of 
graphene plasmons can be combined with the high confinement property of hBN phonon polari- 
tons, allowing for new metamaterials that marry the unique qualities of the two materials. Being 
a hyperbolic material, the hBN provides two completely different regimes of plasmon-phonon 
coupling in the presence of graphene. The spectral dips in spontaneous emission enhancement is 
an example of this coupling where the width of the dips are quite different between the two RS 
bands. Moreover, the contrast for this dip is more than an order of magnitude in both cases with 
the width reaching 20 cm -1 at a frequency of around 800 cm -1 . The width of the spectral dip 
can be controlled by adjusting the Fermi level of graphene and the thickness of the hBN. This can 
provide a powerful route to turning the decay rate on and off via an external knob. In particular, 
this could find application in vibrational spectroscopy 53 and stimulated Raman scattering, with the 
Purcell dip near (Ojo °f either RS band effectively being used as a gate tunable notch filter. Our 
findings might find applications in photodetection, 54 electrically tunable thermal management via 
broadband near field heat transfer 55 and subwavelength imaging. 56,57 The analysis and general 
principles presented in this paper apply equally well to other systems of hyperbolic materials, both 
natural 58 and artificial. 26 
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